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On the Neutrino Flux from Gamma-Ray Bursts 

D. Guetta\ M. Spada\ E. Waxman^ 

ABSTRACT 

Observations imply that gamma-ray bm'sts (GRBs) are produced by the dissipation 
of the kinetic energy of a highly relativistic fireball. Photo-meson interactions of 
protons with 7-rays within the fireball dissipation region are expected to convert a 
significant fraction of fireball energy to > 10^^ eV neutrinos. We present an analysis 
of the internal shock model of GRBs, where production of synchrotron photons and 
photo-meson neutrinos are self-consistently calculated, and show that the fraction of 
fireball energy converted to high energy neutrinos is not sensitive to uncertainties in 
fireball model parameters, such as the expansion Lorentz factor and characteristic 
variability time. This is due in part to the constraints imposed on fireball parameters 
by observed GRB characteristics, and in part to the fact that for parameter values 
for which the photo-meson optical depth is high (implying high proton energy loss to 
pion production) neutrino production is suppressed by pion and muon synchrotron 
losses. The neutrino flux is therefore expected to be correlated mainly with the 
observed 7-ray flux. The time averaged neutrino intensity predicted by the model, 
~ lO^'^'^GeV/cm^ssr, is consistent with the flux predicted by the assumption that 
GRBs are the sources of > 10^^ eV cosmic-rays. 



Subject headings: gamma-rays: bursts-elementary particles-acceleration of particles- 
methods: numerical 



1. INTRODUCTION 

The characteristics of 7-ray bursts (GRBs), bursts of 0.1 MeV-1 MeV photons lasting for a 
few seconds ( [Fishman &: Meegan 1995 ), suggest that the observed radiation is produced by the 



dissipation of the kinetic energy of a relativistically expanding wind, a "fireball," at cosmological 
distance (see, e.g., |Piran 1999 for review). The recent detection of delayed low energy (X-ray to 



radio) emission (afterglow) from GRB sources (see Kulkarni et al. 2000| for review), confirmed 
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both the cosmological origin of the bursts and standard model predictions of afterglows, that 



result from the collision of an expanding fireball with its surrounding medium (see Meszaros 1999 
for review). 

Within the fireball model framework, observed 7-rays are produced by synchrotron emission 
of electrons accelerated to high energy by internal shocks within the expanding wind. In the 
region where electrons are accelerated, protons are also expected to be shock accelerated: Plasma 
parameters in the dissipation region allow proton acceleration to > 10^" eV ( [Waxman 1995a , Vietri 



1995 ; see Waxman 2000| for a recent review). A natural consequence of proton acceleration to 



high energy is the production of a burst of ^ 10^^ eV neutrinos ( [Waxman &: Bahcall 199*7 , 2000), 



produced by the decay of charged pions created in interactions between fireball photons and high 



energy protons. Lower energy neutrinos may be produced by inelastic nuclear collisions (Bahcall 
fc Meszaros 2000| , [Meszaros fc Rees 200q ). 



The characteristic energy of neutrinos produced by photo-meson interactions is determined 
by the relation between the observed photon energy, E^, and the accelerated proton's energy, Ep, 
at the photo-meson threshold of the A-resonance. In the observer frame, 

E^ Ep = 0.2 GeV'^r^, (1) 

where phenomenologically the Lorentz factors of the expanding fireball are F > 10^. For typical 
observed 7-ray energy of 1 MeV, proton energies E'p ~ 2 x 10'' GeV are required to produce 
neutrinos from pion decay, leading to ~ 10^^ eV neutrinos. Both ~ 1 MeV photons and ^ 10^^ eV 
neutrinos are produced in this model during the stage of internal shocks within the expanding 
wind. Much higher energy, ^ 10^^ eV, neutrinos may be produced at a later stage, at the onset of 
fireball interaction with its surrounding medium. Optical-UV photons are radiated by electrons 
accelerated in shocks propagating backward into the fireball plasma. Protons are accelerated to 
high energy in these "reverse" shocks. The combination of low energy photons and high energy 
protons produces ultra- high energy neutrinos via photo-meson interactions, as indicated by Eq. 

The predicted flux of ^ 10^^ eV neutrinos produced in internal shocks is determined by the 
fraction of fireball proton energy lost to pion production. This fraction is determined by the 



number density of photons at the internal dissipation region, and is given by ( Waxman &: Bahcall 



1993,1999) 

UiEp)^0.2minil,Ep/E'p) ^^^'1^ . (2) 

The 7-ray luminosity L^, the photon spectral break energy i?^ (where luminosity per logarithmic 
photon energy interval peaks), the observed variability time At and the wind Lorenz factor F 
are normalized in Eq. (^) to their typical values inferred from observations: = lO^^erg/s, 
E^ = l£^^_MeV MeV, At = lO-^s, and F = lO^-^. The proton break energy, E^, is the threshold 
proton energy for interaction with photons of observed energy E'^, E^ = {2/E^ Mev)-'^2.5 ^ GeV. 
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The rate at which energy is produced as 7-rays by GRBs is similar to the rate of production 
of high energy, > 10^^ eV, protons imphed by the observed ultra-high energy cosmic-ray flux 
([Waxman 1995a| ,b,2000), 



i^,'^ « lO^^erg Mpc-V-^ (3) 



Assuming that GRBs are the sources of observed ultra-high energy cosmic rays, the intensity of 
high energy muon neutrinos and anti-neutrinos implied by Eq. (|2|) is then ( [Waxman fc Bahcall 
1997il999) 

min(l,S^/£;^)GeV/cm2srs. (4) 



E'i,^u ~ 3 X 10"^ 



MEp > El) 



0.2 

Since neutrinos carry ~ 5% of the proton energy, the neutrino break energy is given by 

El^lQ'^l + zr^^^eV. (5) 

-^7,MeV 

We have explicitly introduced here the dependence on source redshift, z, which is due to the fact 
that energies observed on Earth are smaller than those measured at the source redshift by a factor 

The value of /t^, Eq. @, is strongly dependent on F. It has recently been pointed out by 
Halzen and Hooper (1999) that if the Lorentz factor F varies significantly between bursts, with 
burst to burst variations AF/F ~ 1, then the resulting neutrino flux will be dominated by a few 
neutrino bright bursts, and may significantly exceed the flux given by Eq. (^), derived for typical 
burst parameters. This may strongly enhance the detectability of GRB neutrinos by planned 
neutrino telescopes ( [Alvarez- Muniz, Halzen fc Hooper 200q ). 

The main goal of the present work is to determine the allowed range of variation in the 
fraction of fireball energy converted to high energy neutrinos, under the assumption that GRBs 
are produced by internal dissipation shocks in an ultra-relativistic wind. We self-consistently 
calculate the dependence on wind parameters of both 7-ray emission by electrons and neutrino 
production by photo-meson interaction of protons, taking into account synchrotron losses of high 
energy pions and muons and pair-production interaction of high energy photons. We show that 
large burst-to-burst variations in the fraction of fireball energy converted to neutrinos are not 
expected, due to two reasons. First, the observational constraints imposed by 7-ray observations, 
in particular the requirement E"^ > 0.1 MeV, imply that wind model parameters (F, L, At) are 
correlated, and that F is restricted to values in a range much narrower than AF/F ~ 1. Second, 
for wind parameters that imply, through Eq. (|2|), /^r values significantly exceeding 20%, only a 
small fraction of the pions' energy is converted to neutrinos due to pion and muon synchrotron 
losses. 

The model we analyze is described in § |2[ Numerical results are presented and discussed in 
§ ^. Conclusions and implications are summarized in § ^. 
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2. OUTLINE OF THE MODEL 

We consider a compact source, of linear scale Rq ~ 10^ cm, which produces a wind 
characterized by an average luminosity ~ lO^^ergs^^-lO^^ergs"^ and mass loss rate 
M = L^/rj(? (Rq ^ 10^ cm corresponds to three times the Schwarzschild radius of a non-rotating, 
solar mass black hole). At small radius, the wind bulk Lorentz factor, T, grows linearly with 
radius, until most of the wind energy is converted to kinetic energy and T saturates at F ~ 77 ~ 300. 
rj ^ 300 is required to reduce the wind pair-production optical depth for observed high energy, 
> 100 MeV, photons to less than unity. lirj> -q^, ^ {aTLy^/AiTmpC'RoY/'^ = 3 x 10^(L^,53/i?o,6)^'''', 
where = lO^^L^^ssergs"^ and Rq = 10®i?o,6 cm, the wind becomes optically thin at F f« 77=,, < ry, 
and hence acceleration saturates at F ~ 77* and the remaining wind internal energy escapes as 
thermal radiation at ~ 1 MeV temperature. Variability of the source on time scale At, resulting 
in fluctuations in the wind bulk Lorentz factor F on similar time scale, leads to internal shocks 
in the expanding fireball at a radius Ri ~ 2F^cAt = 6 x lO^^F^ 5At_2 cm. If the Lorentz factor 
variability within the wind is significant, internal shocks reconvert a substantial part of the kinetic 
energy to internal energy. It is assumed that this energy is then radiated as 7-rays by synchrotron 
(and inverse-Compton) emission of shock-accelerated electrons. 



The dynamics of our model is described in § |2.1| , and neutrino production calculations are 
described in S 



2.1. Model Dynamics 

In this work, we use an approximate model of the unsteady wind described in the preceding 
paragraph, following Spada, Panaitescu & Meszaros 2000, and Guetta, Spada & Waxman 2000 
(GSWOO). The wind evolution is followed starting at radii larger than the saturation radius, i.e. 
after the shells have already reached their final Lorenz factor following the acceleration phase, and 
the GRB photon flux and spectrum resulting from a series of internal shocks that occur within the 
wind at larger radii are calculated. 

The wind flow is approximated as a set of discrete shells. Each shell is characterized by four 
parameters: ejection time tj, where the subscript j denotes the j-th shell, Lorentz factor Fj, mass 
Mj, and width Aj. Since the wind duration, tw ~ 10 s, is much larger than the dynamical time 
of the source, ~ Rq/c, variability of the wind on a wide range of time scales, < < tyj, is 
possible. For simplicity, we consider a case where the wind variability is characterized by a single 
time scale > t^j > td-, in addition to the dynamical time scale of the source and to the wind 
duration t^. Thus, we consider shells of initial thickness Aj = ctd = Rq, ejected from the source 
at an average rate t~^. 

In GSWOO we have examined the dependence of observed 7-ray flux and spectrum on wind 
model parameters, taking into account both synchrotron and inverse-Compton emission, and 
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the effect of pair production. We have assumed tfiere that the Lorenz factor of a given shell 
is independent of those of preceding shells, and considered various Lorentz factor distributions. 
We have shown that in order to obtain 7-ray flux and spectrum consistent with observations, 
large variance is required in wind Lorentz factor distribution. We therefore restrict the following 
discussion to the bimodal case, where Lorentz factors are drawn from a bimodal distribution, 
Tj = Tm or Tj = Fm ~ ?7* S> T^, with equal probability. The time intervals t^+i — tj are drawn 
randomly from a uniform distribution with an average value of ty. In GSWOO we have considered 
two qualitatively different scenarios for shell masses distribution, shells of either equal mass or 
equal energy, and concluded that observations favor shells of equal mass. Wc therefore restrict the 
following discussion to the equal shell mass case. We note, however, that our conclusions are not 
sensitive to this assumption. 

Once shell parameters are determined, we calculate the radii where collisions occur and 
determine the photon and neutrino emission from each collision. In each collision a forward 
and a reverse shock are formed, propagating into forward and backward shells respectively. 
The plasma parameters behind each shock are determined by the Taub adiabatic, requiring 
continuous energy density and velocity across the contact discontinuity separating the two shells 
(Panaitescu & Meszaros 1999). We assume that a fraction eg (e^) of the protons shock thermal 
energy is converted to electrons (magnetic field). We assume that both electrons and protons are 
accelerated by the shocks to a a power-law distribution, dna/dja c<. "y^^ for particle Lorentz factors 
7a,min < 7a < 7a,max- The maximum Lorenz factor is determined by equating the acceleration 
time, estimated as the Larmor radius divided by c, to the minimum of the dynamical time and the 
synchrotron cooling time. 

The calculation of emitted 7-ray spectrum and flux is carried according to the method 
described in detail in GSWOO, taking into account both synchrotron and inverse Compton 
emission, as well as the effect of the optical thickness due to Thomson scattering on both electrons 
present initially in the flow and those created by pair production. Here we focus on neutrino 
production by photo-meson interactions, following Waxman k, Bahcall 1997. Our method of 
calculation is described below. 



The fraction f-^iEp) of proton energy lost to pion production is estimated as /„■ = min(l, At/t-,^) 
where At is the comoving shell expansion time and is the proton photo-pion energy loss time 
(Waxman & BahcaU 1997), 



2.2. Neutrino production 



1 dEj 
'Ep~~dt 



'p 
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Here (J-k{E) is the cross section for the pion production for a photon with energy E in the proton 
rest frame, ^,{E) is the average fraction of energy lost to the pion, and n[E) is the comoving 
number density of photons per unit energy. In our calculations, we have approximated the integral, 
Eq. (^), by the contribution from the A resonance. This approximation is valid for GRB photon 
spectra (Waxman 2000). 



Photo-meson production of low energy protons, well below E^, requires interaction with high 
energy photons, well above S^. Such photons may be depleted by pair production. For each 
collision we find the photon energy E^ , for which the pair production optical thickness, r^-y, is 
unity. A large fraction of photons of energy exceeding max(ii^^, mgC^) (measured in the shell 
frame) will be converted to pairs, and hence will not be available for photo- meson interaction, 
leading to a suppression of the neutrino flux at low energies. In order to take this effect into 
account, we use /jr = min[l, (At/t7r) min(l, r^^)] for protons interacting with photons which are 
strongly suppressed by pair production. 

Neutrino production is suppressed at high energy, where neutrinos are produced by the 
decay of muons and pions whose lifetime r exceeds the characteristic time for energy loss due to 
synchrotron emission (Waxman & Bahcall 1997, 1999, Rachen & Meszaros 1998). We therefore 
define an effective /tt, fn,eS., as 4/,^ times the fraction of pions' energy converted to muon 
neutrinos. In the absence of pion and muon energy loss, ~ 1/4 of the pions' energy is converted 
to muon neutrinos, since ^1/2 the energy of charged pions is converted to -|- P^. Thus, /7r,efr. 
is the fraction of proton energy that, in the absence of synchrotron losses, leads to approximately 
the same muon neutrino flux as that obtained when synchrotron losses are taken into account. 

In each collision a fraction of the kinetic energy of the colliding shell is converted to a flux of 
photons and neutrinos. The energy that is not lost to photons and neutrinos is converted back to 
kinetic energy by adiabatic expansion of the shell. 



3. RESULTS AND DISCUSSION 

In this section we determine the dependence of /jr^eff.; the fraction of proton energy that in the 
absence of synchrotron losses leads to approximately the same muon neutrino flux as that obtained 
when synchrotron losses are taken into account, on wind model parameters. We adopt electron 
and magnetic field energy fractions close to equipartition, eg = 0.45 and > 0.01, and p = 2. 
The equipartition fractions are required to satisfy eg > 0.1 and eg > 0.01 in order to account for 
observed 7-ray emission (e.g. GSWOO). Moreover, eg values in the range of 0.1 to 0.5 are typically 
derived from GRB afterglow observations (e.g. [Freedman fc Waxman 200 1| ). p ~ 2 is required to 
account for observed 7-ray and afterglow spectra ( [Freedman &: Waxman 2001]) . Figure 1 presents 
contour plots of the photon break energy, E!^, as function of and for = lO^^erg/s and 
and €b = 0.01, = 0.1. The qualitative behavior demonstrated by the contour plot can be 
understood based on the following arguments. The wind energy density at the smallest radii of 
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the internal dissipation shocks is lower for higher values of r^, leading to a decrease of E!^ with 
increasing . This effect sets the upper bound on values. Lowering the value of leads 
to increase in E!^ together with increase in the wind optical depth at the smallest internal shocks 
radii. Once the optical depth exceeds unity, emission is strongly suppressed leading to a rapid 
decrease with Tm of both the fraction of wind energy escaping as radiation and E^ (see GSWOO 
for detailed discussion). 

Figures 2 presents the contour plot for the effective /jr, /7r,efr.) for the case shown in Figure 
1 with eB = 0.01. The values shown are averaged over all internal collisions, weighted by the 
energy converted to high energy protons in each collision. Values of fj^ are shown for several 
neutrino energy ranges, using the approximation Ei, = O.ObEp. The region in Tm-tv plane where 
> 0.1 MeV is bound in Fi gures 2 by the dashed line. The dash-dotted line outlines the region 
in which the fraction of wind energy converted to radiation exceeds 2% (higher fraction is obtained 
for larg er T^rn values). For wind parameters consistent with observed GR.B characteristics, the 
effective value of fj^ at high neutrino energy, E/y > E^, is restricted to values in the range of 
~ 10% to ~ 30%. Figure 3 presents the fraction f^^^'^ of the fireball neutrino flux produced by 
internal shocks which are optically thin to radiation. While at high energy most of the neutrino 
flux is produced in the optically thin regime, where collisions also contribute to the photon flux, 
a large fraction, ~ 50%, of the neutrino flux at low energy is produced by collisions which do not 
contribute to the observed gamma-ray flux. The optically thick collisions provide a significant 
contribution only at low energy, since the large energy density at these collisions lead to strong 
suppression of high energy neutrino production by muon and pion synchrotron losses. 

In figure 4 we present contour plots of /t^, neglecting the effects on neutrino production of 
synchrotron losses and pair-production. approaches unity at low values of F^ for low proton 
(neutrino) energy, and at intermediate values of F^ for high proton (neutrino) energy. The energy 
loss of pions and muons due to synchrotron emission reduces the fraction of pion energy converted 
to neutrinos, and hence suppresses the effective value of /j^, to the range of values shown in Figure 
2. The suppression seen in Figure 4 of /tt at low F^ for high energy protons is due to the fact that 
for low values of F^ most collisions occur at small radii, where synchrotron losses of high energy 



protons prevent acceleration of protons to ultra- high energy ( Waxman 1995a ). Ultra- high energy 
protons are produced in this case only at secondary collisions at large radii, where the photon 
energy density is already low, resulting in smaller values of /t^. 

Figures 5 and 6 present /7r,cfr. for -^w = lO^^erg/s, es = 0.1 and = lO^^erg/s eb = 0.01 
respectively. Comparing to the results presented in Figure 2, we find that the qualitative behavior 
as well as the predicted values of /7r,eff. are not sensitive to wind luminosity and magnetic field 
equipartition fraction (over the range allowed by GRB observations). 

The neutrino break energy inferred from Figures 2, 5 and 6, E^^ ~ 10^^ eV, is somewhat 
higher than estimated by Waxman & Bahcall (1997), -E^ ~ 10^^ eV [see Eq. (^)], since the photon 
energy break in our wind models is typically E!^ ~ 0.1 MeV, rather than ~ 0.5 MeV. Note, that 
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although synchrotron losses suppress neutrino production at high energy (compare Figures 2 and 
4), the increase in 7,^ at high energy (Figure 2) implies that the neutrino flux is flat, with roughly 
equal energy per logarithmic neutrino energy interval above E^, in contrast with the expectation 
of strong suppression of neutrino flux at Ei, ^ 10^^ eV (Waxman & Bahcall 1997, 1999, Rachen 
&; Meszaros 1998). This is due to two effects. First, high energy protons and neutrinos are 
produced in the present model over a very wide range of radii, over part of which synchrotron 
losses are small. Second, the low energy, E^ < E^, photon spectrum obtained in the present 

model, E^{dn-y/dE-y) oc -E^^, is softer than assumed in previous analyses, E^{dn-y/dE^) oc E}^. 
GRB observation suggest that the latter, harder, spectrum is closer to reality, and steepening 
of the spectrum at low photon energy is indeed expected due to effects not included in the 
present analysis, such as inverse-Compton suppression and photospheric emission (see GSWOO for 
discussion). Thus, the present analysis most likely overestimates the value of /7r,efr. for Ey';^ E^. 



4. CONCLUSIONS 

We have studied neutrino production in GRBs, within the framework of the dissipative 
flreball model. We have considered photo-meson production of neutrinos in internal dissipation 
shocks by interaction of accelerated protons with observed 7-rays, taking into account synchrotron 
losses of high energy pions and muons and pair-production interaction of high energy photons. 
Our main results are presented in Figures 2, 5 and 6, showing the dependence on fireball wind 
model parameters of /,r,efF. , the fraction of proton energy that in the absence of synchrotron losses 
leads to approximately the same muon neutrino fiux as that obtained when synchrotron losses 
are taken into account {fj^^cS. is defined as 4/,^ times the fraction of pions' energy converted to 
muon neutrinos). Over the range of wind model parameters, which produce photon break energy 
consistent with observations (bounded by the dashed line in the figures), the value of f-j^^eS. at high 
neutrino energy, Ey > E^, is within the range of 10% to « 30%. The weak dependence of /7r,eff. 
on wind model parameters, in contrast with the strong dependence implied by Eq. (^), is due to 
two reasons. First, for low values of F and At, where large values of /tt are implied by Eq. (^), 
only a small fraction of the pions' energy is converted to neutrinos at high proton energy due to 
pion and muon synchrotron losses (compare figures 2 and 4). Second, the observational constraints 
imposed by 7-ray observations imply that wind model parameters (F, L, At) are correlated. 

Our results imply that GRB neutrino fiux of individual bursts should correlate mainly 
with the bursts' 7-ray flux. An estimate of the time averaged neutrino intensity produced by 
GRBs may be obtained as follows. The average GRB fluence over the energy range observed by 
BATSE, ~ 0.1 MeV to ~ 2 MeV, is 1.2 x 10"^ eV. This corresponds to an average 7-ray intensity 
of E'^i^^ ~ 0.8 X 10^®GeV/ssr cm^. The energy emitted as 7-rays is the energy converted to 
high energy electrons in internal fireball shocks which are optically thin to 7-rays (the electron 
synchrotron cooling time is short compared to the wind dynamical time). Denoting the fraction of 
thermal proton energy converted to electrons in internal shocks by Eg, the expected muon neutrino 
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intensity is 

El^, « i:^il^2i?2$^ « 2-5||-^i^ X lO-^GeV/ssrcm^, (7) 

where /^'^gg- is the value of f-^^cQ, averaged over optically thin collisions only, and f^^^ is the 
fraction of the neutrino flux produced by optically thin collisions (the factor of 2 is due to the 
fact that the electron energy over a logarithmic energy interval is twice that of the photon 
energy per logarithmic photon energy interval). As discussed in §3, the contribution of optically 
thick collisions to neutrino flux is significant only at low energies, see Figure 3. Thus, the ratio 
/w^eff / f^^^ presented in Figure 7 is similar to /7r,cfT. (Figure 2), except at low energy, where it 
exceeds fj^^cs, by a factor of f« 2. 

For eg values implied by observations, typically 0.1 ^ eg ^ 0.5, the neutrino intensity given 
by Eq. (|^, which is based on the observed intensity of GRB 7-ray photons, is similar to that 
given by Eq. (^ derived based on the assumption that GRBs are the sources of ultra-high 
energy cosmic-rays. This similarity reflects the fact that the 7-ray energy generation rate of 
GRBs is similar to the generation rate of high energy protons required to account for the 
observed flux of ultra-high energy cosmic-rays. Note, that we have used many simplifying 
assumptions in our analysis of the variable fireball wind (see § |2|). Hence, the numerical values 
derived by our approximate analysis should not be considered more accurate than the estimates 
given by Eqs. (^^). Nevertheless, the qualitative conclusions regarding the role of synchrotron 
and pair-production suppression and regarding the weak dependence of f-^^eS. on wind model 
parameters are not sensitive to the details of the approximations we have used, and are therefore 
of general validity. We note, in particular, that the value of fj^^cs, does not significantly exceed 
20% also in wind model parameter regions which are outside the parameter region implied by 
observations. 

Finally, we would like to make the following point. Fireball winds with Lorentz factors below 
the range allowed by observations produce low photon luminosity bursts, with characteristic 
photon energies well below 100 keV (see GSWOO). While such bursts may not be detected by 
present photon detectors, and hence we currently have no evidence for their existence, they may 
contribute to the background neutrino intensity. Constraints on the high energy neutrino intensity 
may therefore provide constraints on the existence of such non-GRB fireballs. 
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Fig. 1. — Contour plots of observed photon spectral break energy (where photon luminosity per 
logarithmic photon energy interval peaks), logio(-E'7[eV]) as function of wind variability time ty and 
minimum Lorentz factor F^, for = lO^^erg/s and 63 = 0.01 (left panel), es = 0.1 (right panel). 
The wind maximum Lorentz factor is Fm = 2500, and the source is assumed to lie at z = 1.5. 
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Fig. 2. — Contour plots of fn,cS., the effective value of /t^ (defined as 4/,^ times the fraction of 
pions' energy converted to muon neutrinos), as function of wind variability time and minimum 
Lorentz factor T^n, for wind luminosity = lO^^erg/s and cb = 0.01. The four panels 
correspond to four observed neutrino energy bins, clockwise from top left: lO^^eV < E^, < lO^^eV, 
lO^^eV < < IQi^eV, lO^^eV < E^ < lO^^eV, lO^^eV < E^ < lO^^eV. We have used the 
approximate relation Ey = 0.05ii^p between neutrino and proton energy, and assumed a source 
redshift z = 1.5. The region in Pi^^^^t,; plane where E^, ^0.1 AlcV is bound by the dashed lines. 
The dash-dotted lines outline the region in which the fraction of wind energy converted to radiation 
exceeds 2% (higher fraction is obtained at larger values). 
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Fig. 3. — Contour plots of /* , the fraction of the neutrino flux produced by optically thin 
collisions, as function of wind variability time ty and minimum Lorentz factor T^, for the case 
shown in Figure 2. 
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Fig. 4. — Contour plots of /tt, obtained when the effects on neutrino production of synchrotron 
losses and pair-production are neglected, for the case shown in Figure 2. 
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Fig. 7. — Contour plots of the ratio f^^g / f^}^^^ , as function of wind variability time ty and minimum 
Lorentz factor F^, for the case shown in Figure 2. f^^^ is the value of /,T,efr averaged over optically 
thin collisions only and Z^'^™, the fraction of the neutrino flux produced by optically thin collisions, 
is shown in Fig. 3. 



